Evidence for Spinodal Decomposition in Nuclear Multifragmentation 
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Multifragmentation of a "fused system" was observed for 
central collisions between 32 MeV/nucleon 129 Xe and nat Sn. 
Most of the resulting charged products were well identified 
thanks to the high performances of the INDRA 4-zr array. 
Experimental higher-order charge correlations for fragments 
show a weak but non ambiguous enhancement of events with 
nearly equal-sized fragments. Supported by dynamical cal- 
culations in which spinodal decomposition is simulated, this 
observed enhancement is interpreted as a "fossil" signal of 
spinodal instabilities in finite nuclear systems. 
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The decay of highly excited nuclear systems through 
multifragmentation is, at present time, a subject of great 
interest in nucleus- nucleus collisions, and many efforts 
are made to understand the underlying physics. While 
this process has been observed for many years, its exper- 
imental knowledge was strongly improved only recently 
with the advent of powerful 47r devices which authorize 
careful selections of well defined fused systems undergo- 
ing multifragmentation. Many theories have been de- 
veloped to explain multifragmentation (see for example 
ref. jj| for a general review of models). One theory, in 
particular, is the concept of multifragmentation which 
considers volume instabilities of the spinodal type. In- 
deed, during a collision, a wide zone of the nuclear matter 
phase diagram may be explored and the nuclear system 
may enter the liquid-gas phase coexistence region (at low 
density) and even more precisely the unstable spinodal 
region (domain of negative incompressibility) J2). Thus, 
a possible origin of multifragmentation may be found 
through the growth of density fluctuations in this un- 
stable region. Within this theoretical scenario a breakup 



into nearly equal-sized "primitive" fragments should be 
favored in relation with the wave-lengths of the most un- 
stable modes present in the spinodal region ||. How- 
ever this simple picture is expected to be strongly blurred 
by several effects: the beating of different modes, even- 
tual coalescence of nascent fragments, secondary decay 
of excited fragments and mainly the finite size of the sys- 
tem Pj. Therefore only a weak proportion of multifrag- 
mentation events with nearly equal-sized fragments is ex- 
pected. Note that surface instabilities developed from pe- 
culiar multifragmenting geometrical structures (bubbles, 
disk-like), can also produce such events In this letter 
we investigate, using a very sensitive method of charge 
correlations, the occurrence of equal-sized fragment par- 
titions in a selected sample of experimental events, cor- 
responding to a fused system undergoing multifragmen- 
tation. A comparison is made with predictions of 3D 
stochastic mean-field simulations of collisions which take 
into account the dynamics of the most unstable modes 
in the spinodal region. 

We have experimentally studied 

the reaction 129 Ae + nat Sn at 32 MeV/nucleon. The 
measurements were performed at the GANIL accelerator 
using the 47r multidetector for charged reaction products 
INDRA H . Detailed information on the experiment and 
on the selection of fused events can be found in Ref. . 
Let us just recall that the selection was performed by re- 
quiring the detection of a significant fraction (> 80%) of 
the charge of the system; these selected events are called 
complete events. Reaction products with charge Z > 5 
were defined as fragments. Finally the preferred direc- 
tion of emission of matter in the center of mass of the 
reaction (flow angle) was determined from the calcula- 
tion of the energy tensor of fragments, and the requirc- 
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ment was made that this angle be larger than 60° 
The main argument underlying the chosen selection is 
that while a fused system should be present at all flow 
angles, binary dissipative collisions should vanish when 
this angle is large, giving way to an almost pure "fusion" 
phenomenon Q.The present selection corresponds to ^25 
mb, and the total cross-section for multifragmentation of 
compact shape systems is estimated to 100 mb. Figure |l| 
exhibits how complete events populate the flow angle do- 
main as a function of their measured total kinetic energy 
(i.e. cm. kinetic energy of all detected fragments and 
light charged particles). 
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FIG. 1. Wilczynski diagram for complete events: correla- 
tion between total measured c.m. kinetic energy and flow 
angle 9fi a w 

An unambiguous charge identification is essential for 
building unbiased charge correlation functions: in this 
experiment all fragments with charge Z<20 are identi- 
fied with a resolution of one charge unit. The correlation 
method used, called "higher order correlations" was pro- 
posed a few years ago in ref |jTo|| . It is appropriate to 
search for weak signals and its originality consists of the 
fact that all information on fragments of one event is 
condensed in two variables to construct the charge cor- 
relation. This quantity is defined by the expression: 

Y(AZ,<Z>) 

Y'(AZ, < Z >) M { > 

Here, Y(AZ, < Z >) is the yield of selected events with 

< Z > and AZ values; M is the fragment multiplicity, 

< Z > denotes the average fragment charge of the event: 

1 M 

< Z >=MZZ Z > ( 2 ) 

i=i 

and AZ is the standard deviation, defined by: 



Az = ^Mn$>«- < z » 2 ( 3 ) 

The denominator Y'(AZ, < Z >) which represents the 
uncorrelated yield is built, for each fragment multiplicity, 
by taking fragments in different events of the selected 
sample. The number of uncorrelated events is chosen 
large enough (10 3 per true event) to strongly reduce their 
contribution to statistical error. With such a correlation 
method, if events with nearly equal-sized fragments are 
produced, we expect to see peaks appearing in the first 
AZ bin (0-1). 




FIG. 2. Experimental higher-order charge correlations for 
selected fusion events in 32 MeV/nucleon Xe+Sn collisions. 

Higher order correlation functions for selected experi- 
mental events are shown in fig. |^. We observe peaks in 
the bin AZ = 0-1 for each fragment multiplicity. From 
M=3 to M=6 the maximum is moving from < Z >=19 
to < Z >=10. When multiplied by the multiplicity these 
lead to a roughly constant value close to the most prob- 
able total fragment charge for all fused events. It was 
moreover verified that the kinetic energy of fragments in 
these events does not differ (within statistical error) from 
average kinetic energies of fragments for other selected 
events. 

We can now estimate whether the enhancement of 
events with equal-sized fragments is statistically signif- 
icant and quantify their occurrence. In this aim we built 
charge correlations for all events, whatever their mul- 
tiplicity, by replacing the variable < Z > by Z to t — 
Mx < Z >. For this compact presentation uncorre- 
lated events are built and weighted in proportion to real 
events of each multiplicity. For each bin in Z to t, fixed 
at six atomic number units, an exponential evolution of 
the correlation function is observed from AZ=7-8 down 
to AZ=2-3. This exponential evolution is thus taken as 
"background" to extrapolate down to the first AZ bin. 
Higher order correlation functions for the first bin in AZ 
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are displayed in Fig. || with their statistical errors; the 
solid line corresponds to the extrapolated "background" . 
All events corresponding to the three points whose error 
bar is fully located above this line correspond to a statis- 
tically significant enhancement of equal-sized fragment 
partitions. The probabilities that these values which 
are higher than the background simply arise from sta- 
tistical fluctuations are 0.048 (Z tot =45), 0.038 (Z tot =51) 
and 0.022 (Z tot =57). The number of significant events 
amounts to 0.1% of selected fusion events. 




FIG. 3. Higher-order charge correlations: quantitative re- 
sults for experimental data (upper panel) and simulations 
(lower panel). Symbols indicate the events where AZ = 0-1, 
curves show the background defined extrapolating AZ > 2 
(see text). Vertical bars correspond to statistical errors as- 
suming independent measurements and horizontal bars define 
Ztot bins. 



C 4.S 
O 4 

a 3 

£2.5 

o 24 
Qi.s 

M 
0.5 h 




M=3 



3^ 



M=4 





FIG. 4. Experimental higher-order charge correlations for 
very dissipative "binary" collisions in 32 MeV/nucleon Xe+Sn 
collisions (see text). 

To test if these equal-sized fragment partitions were 
specific - or not - of the selected fused events, we have 
studied a data sample dominated by very dissipative col- 
lisions of binary type. This is done by selecting in the 
Wilczynsky diagram for complete events (see fig. |l|) to- 
tal cm. energy lower than 1200 MeV with a flow angle 
0fiow< 30°. Within this new selection we estimate that 
about 9xl0 3 fused events are mixed with about 3.6xl0 5 
binary (or eventually strongly deformed, or both) events. 
The correlation functions are presented in fig. No sta- 
tistically significant signals are observed. The small peak 
for M=6 corresponds to 0.0008% of selected events. This 
is an indication that only fused events are characterized 
by an enhanced production of nearly equal-sized frag- 
ments. 

In ref. [||, it was experimentally observed that charge 
distributions for fusion/multifragmentation events were 
independent of the total mass of the incident nuclei (248 
or 393 u), while fragment multiplicities scaled as their 
total charge. These properties were the first hint of a 
bulk effect for producing fragments. 

Dynamical stochastic mean-field simulations j[l 12 
were performed for head-on collisions only, thus neglect- 
ing shape effects. Spinodal decomposition is simulated 
using the Brownian one-Body (BoB) dynamics [fi"3| p"5| , 
which consists in employing a Brownian force in the ki- 
netic equations. The ingredients of the simulations are 
as follows. The self-consistent mean field potential fil| 
chosen gives a soft equation of state {Koo— 200 MeV) and 
the finite range of the nuclear interaction is taken into ac- 
count using a convolution with a gaussian function with 
a width of 0.9 fm [[l7|. The addition of a term propor- 
tional to Ap in the mean-field potential allows to well- 
reproduce the surface energy of ground-state nuclei [|l8| . 
This is essential in order to correctly describe the ex- 
pansion dynamics of the fused system. In the collision 
term a constant a nn value of 41 mb, without in-medium, 
energy, isospin or angle dependence is used filj|| . In the 
following step the spatial configuration of the primary 
fragments, with their excitation energies as produced by 
BoB, was taken as input in the SIMON code (2^] to 
follow the fragment deexcitation while preserving space- 
time correlations. Finally the events were filtered to ac- 
count for the experimental device. These complete sim- 
ulations well reproduce multiplicity and charge distribu- 
tions of fragments and their average kinetic energies jll| . 
To refine the comparison, higher-order charge correla- 
tions are calculated for the simulated events. Although 
all events in the simulation arise from spinodal decom- 
position, only a very small fraction, as suspected for fi- 
nite systems [|| , exhibits final partitions with equal-sized 
fragments. This fraction is not strongly modified by the 
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de-excitation of fragments; as an example Z=15 primary 
fragments produced in BoB simulations (with their mass 
and excitation energy distributions) lead to a secondary 
Z distribution centered at Z=14 with a standard devia- 
tion of 0.6. Finally the proportion of statistically signif- 
icant equal-sized fragment partitions is similar (0.15%) 
to the experimental one. A detailed quantitative com- 
parison is displayed in fig. ^|. The similarities between 
experimental and calculated events allow to attribute all 
fusion-multifragmentation events to spinodal decomposi- 
tion. The peaks observed near AZ=0 in the higher-order 
charge correlations are thus fossil fingerprints of the par- 
titions expected from spinodal decomposition. It was 
previously observed that multiplicity and Z distributions, 
and the average kinetic energy of fragments were also 
well reproduced with the statistical model SMM |^l],^H . 
It is interesting to underline that in this sample of SMM 
events no partition with AZ <2 exists. Note also that 
fragment velocity correlations |23f| are perfectly repro- 
duced either in BoB simulations, or in the breakup of 
systems with elongated shapes. As none of these config- 
urations is strongly necked or close to a bubble, we can 
exclude surface instabilities as a possible source of equal- 
sized fragments, in favor of bulk instabilities flu] f24| , |25| ] . 

In conclusion, we have investigated charge correlation 
functions for fusion events which undergo multifragmen- 
tation. We have found a weak but unambiguously en- 
hanced production of events with equal-sized fragments. 
Supported by theoretical simulations we interpret this 
enhancement as a signature of spinodal instabilities as 
the origin of multifragmentation in the Fermi energy do- 
main. Spinodal instabilities are thus shown for the first 
time in a finite system. Moreover the occurrence of spin- 
odal decomposition shows the presence of a liquid-gas 
coexistence region and gives a strong argument in favor 
of the existence of a first order liquid-gas phase transition 
in finite nuclear systems. 
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